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ABSTRACT 



Aims. We gather optical spectra of 8 long-duration GRB host galaxies selected from the archival data of VLT/FORS2. We investigate 
whether or not Wolf-Rayet (WR) stars can be detected in these GRB host galaxies. We also estimate the physical properties of GRB 
host galaxies, such as metallicity. 

Methods. We identify the WR features in these spectra by fitting the WR bumps and WR emission lines in blue and red bumps. We 
identify the subtypes of the WR stars, estimate the numbers of stars in each subtype, and calculate the WR/O star ratios. The (O/H) 
abundances of GRB hosts are inferred from both the electron temperature (T e ) and the metallicity-sensitive strong-line ratio {Rm), for 
which we break the R23 degeneracy. We compare the environments of long-duration GRB host galaxies with those of other galaxies 
in terms of their luminosity (stellar mass)-metallicity relations (L - Z, M, - Z). 

Results. We detect WR stars in 5 GRB host galaxies with spectra of relatively high signal-to-noise ratios (5 /N). In the comparison 
of L — Z, M, — Z relations, we show that GRB hosts have lower metallicities than other samples of comparable luminosity and stellar 
mass. The presence of WR stars and the observed high WR/O star ratio, together with the low metallicity, support the "core-collapsar" 
model and imply that we are witnessing the first stage of star formation in the host regions of GRBs. 

Key words. Gamma rays: bursts - Stars: Wolf-Rayet - Galaxies: abundances - Galaxies: fundamental parameters 



1. Introduction 

Gamma-ray bursts (GRBs), the most energetic events 
in the Universe, were discovered a ccidentally in the 
1960s dKlebesadel. Strong & Olsonl Il973l) . Since then, GRBs 
have been the targets of intense research. However, the mech- 
anism behind these bursts and the identity of their progenitors 
remain disputed. 

The durations and spectral properties of GRBs suggest 
a classification of short (of duration < 2 s and hard spec- 
trum) and long bursts (of duration > 2 s and soft spectrum) 
dKouveliotou et al.|[l993t lHartmannll2005l) . Short GRBs are be- 
lieved to originate from the merger of compact binaries such 
as double neutron star binaries (NS-NS) and black hole - 
neutron star binaries (BH-NS). For long GRBs, the favored 
"core-collapsar" model starts from the idea of a rapidly rotat- 
ing, massive star that has undergone extre me gravitational col- 
lapse and formed a central black hole ( Wooslev et al.l 119931: 
iMacFadven & WooslevlfT^l iKlose et al]T2004l) . According to 
the collapsar model, the prompt energetic structure of a long 
GRB is the result of energy dissipation by internal, relativis- 
tic shocks, which may last seconds or minutes, at a radius of 
about 10 14 cm from the center of the collapsed star (Hartmann 
I2005L and references therein). Moreover, the association be- 
tween GRBs and supernovae (SNe) indicates that in many cases 
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the parent SN population of GRBs is formed by peculiar type 
IbcSNe. 

Wolf-Rayet (WR) stars are naturally considered to be the 
most favored candidates of lo ng duration G RB progenitors. In 
the stellar evolution model of Hirschi et al.l d2005l) . WR stars, 
which are massive short-lived stars, satisfy the main criteria for 
GRB production of black hole formation, loss of hydrogen-rich 
envelope, and sufficient angular momentum to form an accre- 
tion disk around the black hole. According to this model, a 
lower limit to the metallicity of subsolar value (typically be- 
tween Zsmc and Zlmc, i-e-, Z ~ 0.2-0.4 Z ) is also a criterion 
for GRB production. However, magn etic-field breaking poses 
some difficulties in producing GRBs (Petr ovic et al .1120051) . To 
solve this problem, one possible scenario involves assuming that 
the star of lower metallicity rotates so rapidly that mixing oc- 
curs and the star chemically evolves homogeneously without a 
hydrogen envelope. Moreover, its lower metallicity (typically 
Z < 0.05 Zq) causes low mass loss and therefore the retention 
of a high angula r momentum. Those conditions are necessary t o 
produce a GRB dWooslev & Hegerl2006l:lYoon & Langej2 005). 
Therefore, it is imperative to confirm the presence of WR stars 
and determine more accurately metallicities in the region of 
GRBs. 

The subtype and number of WR stars, and the relative WR/O 
star number ratio, are related to the star-forming activity and the 
starburst duration in galaxies. We can achieve a deeper under- 
standing of the evolutionary paths of long-duration GRB progen- 
itors by detecting the WR populations within GRB host galax- 
ies. Moreover, WR evolutionary models suggest that metallicity, 
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one of the important diagnostics of the evolut ionary histories of 
galax ie s, affects the properties of W R stars dSchaerer & Vaccal 
119981) . ICrowther & Hadfieldl d2006l) investigated the effect of 
metallicity on the WR/O star ratio. They found that the WR/O 
star ratio decreases with metallicity. However, high WR/O star 
ratios are found in host galaxies of both GRB 980425 and GRB 
020903, which have low metallicities (Ham mer et al.l 12006). 
The WR/O star ratio - metallicity relation for other GRB hosts is 
still unknown. 

The first detection of a counterpart to a GRB at optical and 
X-ray was achieved on 28 February 1997 (van Paradiis et al. 
|1997|) . Since the discovery of GRB afterglows, GRB host galax- 
ies and their redshifts have been identifiable. Between 1997 and 
2007, 588 GRBs were detected. Among them, 325 GRBs had X- 
ray afterglows, 220 GRBs had optical afterglows, and 55 GRBs 
had radio afterglows Q. 

The study of GRB host galaxies is very important to the un- 
derstanding of the physical properties of GRB regions and the 
nature of GRB progenitors. Although a large number of host 
galaxies can be identified, most of them are too faint to be ob- 
served even using the largest telescopes in the world. So far, only 
a few dozens of host galaxies have been spectroscopically ob- 
served (e.g., GHostS). Moreover, the number of GRB hosts that 
have been intensively studied spectroscopically is even less. 

Evidence from photometric and spectroscopic observa- 
tions shows that the host galaxies of long-duration GRBs 
are mostly faint, blue, low-mass, star - forming galaxies with 
low metallicities dSokolov et al l I200U iLe Floc'h et all [2003: 
Fvnbo et al.l 120031; ICourty et al.l 120041 iProchaska et al.l 12 004 ; 
Christens en et al.1 |2004t ICharvetaU l2002 | : | G orosabe l et al 
20051: iFvnbo et all 120061 Iwiersema et al.l 120071: iKewlev et al 



Table 1. Basic information of the GRB host galaxies 



2007E iLevesque et all 120091: ISavaglio et al. 2009, hereafter 
Savaglio09). In contrast, the host galaxies of short-duration 
GRBs mostly have higher lum inosities and higher metallicities 
than long-duration GRB hosts (Ber gerll2008l) . 

We study a sample of 8 long-duration GRB hosts with high 
quality spectroscopic observations. Firstly, we try to detect the 
WR features in these GRB hosts then study the physical proper- 
ties of GRB host galaxies, such as their metallicities and lumi- 
nosity (stellar mass)-metallicity relations (L-Z, M*-Z). This pa- 
per is organized as follows. The sample selection and flux mea- 
surements are performed in Sect. [2] In Sect. [3] we describe the 
identification of the WR features in the spectra of GRB hosts. 
In Sect. |U the physical properties of GRB host galaxies are dis- 
cussed. The discussion and conclusions are presented in Sect. 
Throughout the paper, we adopt the ACDM cosmological model 
(H = 70 km s -1 Mpc~', Q M = 0.3, and Q A =0-7), and the initial 
mass function (IMF) proposed by Sal peterl (Il955l) . All compar- 
isons performed in this work include a normalization to the same 
ACDM cosmological model and IMF. All magnitudes in this pa- 
per are in the Vega system. 



2. Data reduction and measurements 

2.1. The sample selection and data reduction 

One of the main goals of this work is to find evidence of WR 
stars in GRB hosts by means of optical spectroscopic analysis. 
We searched archival data of VLT/FORS4J to obtain the spectra 
of GRB host galaxies at z < 1, which ensures that a sufficient 



GRB 




RA a 




DEC 




z 


E G (B - V) b 


Type 


GRB 980703 


23 


: 59 : 


07 


08 : 


33 : 


36 


0.966 


0.061 


long 


GRB 990712 


22 


: 31 : 


50 


-73 


: 24 


: 29 


0.433 


0.033 


long 


GRB 020405 


13 


: 58 : 


03 


-31 


: 22 


: 22 


0.691 


0.055 


long 


GRB 020903 


22 


: 49 : 


25 


-20 


: 53 


: 59 


0.251 


0.033 


long 


GRB 030329 


10 


: 44 : 


50 


21 : 


31 : 


17 


0.168 


0.025 


long 


GRB 031203 


08 


: 02 : 


28 


-39 


: 51 


: 04 


0.105 


1.040 


long 


GRB 060218 


03 


: 21 : 


39 


16 : 


52 : 


02 


0.034 


0.140 


long 


GRB 060505 


22 


: 07 : 


01 


-27 


:48 


: 56 


0.089 


0.021 


long? 



Notes : (a) Coordinates at the 2000.0 epoch, (b) Galactic extinction 
dSchlegel et alfl 998) 



number of emission lines are included in the spectral coverage. 
These spectra taken with the 600B, 600RI, 600Z grisms (R ~ 
1300), and 300V grism (R ~ 400) of FORS2 were adopted, be- 
cause they represent a good compromise between relatively high 
spectral resolution, which allows to resolve the WR features, and 
relatively high sensitivity, which provides spectra of high qual- 
ity. 

Data reduction and extraction were performed using a set of 
IRAF3 procedures developed by our team, which can simultane- 
ously reconstruct the spectra and the sky counts of the objects. 
For some GRBs associated with supernovae, the early-time spec- 
tra were carefully checked to avoid the contamination of super- 
novae, which increases the difficulty in identifying WR features 
and the error in the emission line flux. Spectra containing con- 
tamination by supernovae were then removed. 

Wavelength calibration was performed using the HeNeAr 
lamp spectra. Flux calibration was achieved by observing 
spectrophotometric stars with the same grism and the same 
slit width. Galact i c exti nction was adopted from NEE0 and 
ISchaerer & Vaccal l[l998). All spectra were transformed into the 
rest-frame. The spectra taken with different grism, but the same 
resolution and the same slit width were merged to widen the 
wavelength coverages (e.g., GRB 980703; GRB 020903). When 
merging spectra, the exposure times were applied as weighting. 
The sample has 8 objects including 7 long-duration GRB hosts 
and 1 possible long-duration GRB host (GRB 060505). The host 
of GRB 060505 was spatially resolved by VLT/FORS2 spec- 
troscopic observations, therefore the GRB region and the en- 
tire host galaxy were studied, separately (see lThone et al.l f2008. 
for the imaging). The detailed discussion about GRB 060505 is 
given in Sect. [5] 

Table[T]lists those 8 GRB hosts in our sample with the name, 
coordinate at 2000 epoch, redshift, and burst type. The detailed 
observational information (date, exposure time, seeing, grism, 
and program) for the objects is given in Table [2] 

2.2. Flux measurements of emission lines 

To measure accurate emission line fluxes, the continuum and 
absorption lines should be subtracted from the spectrum care- 



For the most complete list of GRBs, see the URL: 
http://www.mpe.mpg.de/~jcg/grbgen.html maintained by J. Greiner. 
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U 1 Date 




Exposure time(second) 


seeing a (") Grism 


Pro 


'ram ID 


GRB980703 














2004.07.15 - 


17 


8 x 1200 


0.75 


600Z 


073.B 


- 0482(A) 


2004.07.16 - 


17 


8 x 1200 


0.75 


600RI 


073.B 


- 0482(A) 


GRB990712 














2005.07.05 - 


06 


4 x 1800 


0.59 


300V 


075.D 


-0771(A) 


GRB020405 














2004.07.15 - 


16 


9 x 1200 


0.70 


600RI 


073.B 


- 0482(A) 


GRB020903 














2004.07.15 




6 x 1200 


0.55 


600B 


073. B 


- 0482(A) 


2004.07.16 




6 x 1200 


0.83 


600RI 


073.B 


- 0482(A) 


GRB030329 














2003.04.10 




6x600 


1.18 


300V 


071.D 


- 0355(B) 


2003.04.17 




2x300 


0.71 


300V 


071.D 


-0355(B) 


2003.04.22 




1 x300 


0.71 


300V 


071.D 


- 0355(B) 


2003.05.01 




2x600 


0.51 


300V 


071.D 


- 0355(B) 


2003.06.19 




3 x900 


0.75 


300V 


271.D 


- 5006(A) 


GRB031203 














2003.12.20 




2 x 2700 


0.50 


300V 


072.D 


- 0480(A) 


2004.09.20 




2 x 1800 


0.60 


300V 


073.D 


- 0255(A) 


2005.03.28 




4 x 1200 


0.69 


300V 


073.D 


-0418(A) 


2005.04.09 




3 x 1200 


1.26 


300V 


075.D 


-0771(A) 


GRB060218 














2006.11.27 




1 x 3100 


0.97 


300V 


078.D 


- 0246(A) 


2006.12.19 - 


20 


3x3150 


0.72 


300V 


078.D 


- 0246(A) 


GRB060505 














2006.05.23 




2x 1800+ 1327 + 600 


0.81 


300V 


077.D 


-0661(B) 



Notes: (a) Seeing is the average value in the observing night. 



fully. To do t his, the STARLIGHT0, a sp ectral synthesis code, 
developed by ICid Fernandes et all ( 20051) . is used to fit the ob- 



served spectrum. STARLIGHT can model the continuum and 
stellar absorption lines using a linear combination of N* simple 
stellar population (SSP) templates from the evoluti onary pop- 
ulation synthesis code of lBruzual & Chariot I (120031) . In the fit- 
ting, we use a base of 45 SSPs, which includes 15 ages (0.001, 
0.003, 0.005, 0.01, 0.025, 0.04, 0.1, 0.28, 0.64, 0.9, 1.4, 2.5, 5, 
10, 13 Gyr) at 3 metallicities (0.2, 1, 2.5 Z Q ). The Galactic ex- 
tinction law of ICardelli et all d 1989b (CCM) was adopted in the 
STARLIGHT fitting. An example of spectral fitting is shown in 

Emission line fluxes were measured manually using SPLOT 
task in IRAF. The errors originated mainly from three sources: 
the first is the uncertainty in fitting both the continuum and the 
stellar absorption line; the second is the uncertainty in the flux 
measurement; the third is the Poisson noises from both sky and 
objects, which dominate the error budgets. Emission line fluxes 
corrected for Galactic extinction are given in Table [3] 



2.3. Dust extinction 

The dust extinction can be estimated from the Balmer-line ratios 
(Ha/HjS, HyfHB, and Ha/Hy). In our dust extinction calculation, 
we adopted Case B recombination with an electron density of 
100 ctrr 3 and a temperature of 10 000 K. The predicted intrin- 
sic ratio is 2.86 fo r I (Ua)/I (HB), and 0.466 for / (H r )// (H y 8) 
(lOsterbrockl 1 1 989b . We applied the Galactic extinction law of 
CCM for R v = 3.1 where R v is the ratio of A v to E(B - V) 




5000 6000 
Wavelength (A) 



Fig. 1. The spectrum of host galaxy of GRB 031203. The spec- 
tral synthesis fitting is plotted in red line. The strong emission 
lines are marked. 



http://www.starlight,uf sc.br 



dSeatonll979l) . The values of dust extinction are listed in Table|4] 
The 2 extinction values for each object derived from the ra- 
tios of Ha/HB and HyfHB are consistent with each other. The 
E(B - V) derived from Ha/HB was adopted in this study when 
available, because the higher S /N ratios of these two lines can 
ensure smaller uncertainties in the results. The E(B - V) derived 
from HyfHB was adopted for GRB 980703 and 020405, because 
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Table 3. Emission line fluxes of GRB host galaxies (in units of 1CT 17 erg s _I cirT 2 , corrected for Galactic extinction) 



Ion 


980703 


990712 


020405 


020903 


030329 


031203 


060218 


060505" 


060505" 



[OII]^3727 


24.28 ±0.11 


23.26 ±0.15 


11.53 ±0.16 


11.48 ±0.36 


19.62 ± 0.32 


943.07 ± 4.80 


224.25 ±1.01 


20.96 ± 0.82 


50.33 


± 


1 


88 


HS 


1.66 ±0.22 


1.95 ±0.09 


- 


1.76 ±0.24 


- 


204.83 ±4.41 


20.11 ±0.57 


- 










Hy 


3.1 ±0.19 


4.9 ±0.10 


1.9 ±0.09 


3.94 ± 0.29 


4.41 ±0.10 


499.73 ± 1.10 


42.02 ± 0.67 


3.79 ± 0.30 










[OIII]/14363 




0.92 ± 0.08 




0.55 ± 0.07 


0.71 ± 0.05 


79.10 ± 1.25 


4.77 ± 0.30 












[FeIII]/!4658 


0.16 ±0.09 
























HeIU4686 






0.05 ± 0.03 


0.39 ± 0.06 




1.20 ±0.92 


0.81 ±0.31 












[ArIV]i4711 






0.07 ±0.10 


0.80 ± 0.25 




4.90 ±3. 11 


2.12 ±0.84 












[ArIV]/i4740 


1.01 ±0.05 






0.41 ±0.10 




3.81 ±2.91 


2.28 ±0.91 












w 


6.79 ± 0.30 


11.29 ±0.07 


5.36 ±0.10 


8.58 ± 0.29 


9.80 ±0.16 


1135.37 ±5.21 


91.43 ±0.65 


8.90 ± 0.32 


21.02 


± 





51 


[OIII]i4959 


4.51 ±0.41 


17.37 ±0.11 


5.19 ±0.06 


15.56 ±0.31 


11.74 ±0.20 


2404.78 ± 5.32 


108.4 ± 0.48 


8.47 ± 0.24 


12.47 


+ 





32 


[OIII]i5007 


14.39 ±0.33 


47.28 ± 0.09 


18.33 ±0.13 


44.11 ±0.33 


30.53 ± 0.26 


7238.76 ± 4.68 


291.23 ±0.66 


20.48 ± 0.48 


27.80 


+ 


2 


30 


Ha 




40.34 ± 0.32 




26.03 ± 0.21 


30.94 ± 0.27 


3636.63 ± 4.96 


261.51 ±0.54 


31.49 ± 0.44 


114.8 


± 


3 


66 


[NII]/16583 




< 0.41 ± 0.20 




0.77 ±0.10 


0.29 ±0.15 


109.04 ±4.21 


10.73 ±0.32 


1.79 ±0.09 


29.62 


± 





29 


[SIIJ/16716 








2.14 ±0.17 


3.82 ± 0.09 


129.55 ±4.33 


15.4 ±0.36 


3.46 ±0.13 


26.36 


+ 





51 


[SII]^6731 








1.28 ±0.09 


1.89 ±0.08 


94.63 ± 2.89 


11.56 ±0.41 


2.80 ± 0.20 


13.98 


+ 





19 



Notes: (a) GRB site, (b) The entire host galaxy. 



Table 4. Dust extinction of GRB host galaxies 



the subtypes and numbers of WR stars, and finally summarize 
the results of these GRB host galaxies. 



GRB 

980703 

990712 

020405 

020903 

030329 

031203 

060218 

060505° 

060505* 



Ha/HyS Hy/H/J A v (Ha/Hfi) A v (Hy/HS) E HG (B-V) 



3.57 ±0.41 

3.03 ±0.91 
3.16 ±0.57 
3.20 ± 1.79 
2.86 ± 1.64 
3.54 ±3.59 
5.46 ±4.51 



0.46 ± 0.23 
0.43 ±0.11 
0.35 ±0.10 
0.46 ± 0.32 
0.45 ±0.12 
0.44 ± 0.25 
0.46 ± 0.34 
0.43 ± 0.72 



0.55 ± 0.03 

0.15 ±0.09 
0.24 ± 0.05 
0.28 ± 0.04 
0.00 ± 0.02 
0.53 ± 0.08 
1.60 ±0.31 



0.14 ±0.53 
0.50 ±0.15 
1.93 ±0.36 
0.10 ±0.57 
0.25 ± 0.20 
0.40 ±0.12 
0.10 ±0.12 
0.63 ± 0.35 



0.05 ± 
0.18 ± 
0.62 ± 
0.05 ± 
0.08 ± 
0.09 ± 
0.00 ± 
0.17 ± 
0.52 ± 



0.17 
0.01 
0.12 
0.03 
0.02 
0.01 
0.01 
0.03 
0.10 



Notes: (a) GRB site, (b) The entire host galaxy. 



Ha is not available in their spectra. We then corrected emission 
line fluxes for dust extinction. 



3. Wolf-Rayet bump identification 

The main WR features often seen in optical spectra of galax- 
ies are two characteristic broad emission line clusters. One is 
in a blue part o f the spectrum a t around 4600-4680 A (here- 
after blue bump.lAllen et al.ll976t iKunfh & Jouberl 19851: IContil 



1991; Schaerer&Vaccal 119981) . The other is in a red region 
around 5650-5800 A (herea fter red bump. lKunfh & Schildl 1986t 
iDinerstein & Shieldslll986l) . The blue bump is actually a blend 
of some broad WR lines, such as N V 4605, 4620 A, N III 
4634, 4640 A, C III/C IV 4650, 4658 A, and He II 4686 A, and 
some nebular emission lines superimposed on the bump, such 
as [Fe III] 465 8 A, He II 4686 A He I + [Ar IV1 47 1 1 A, and 
[Ar IV] 4740 (IGuseva et all 120001: 1 Izotov et al.l[l998l The red 
bump is usually much weaker than the blue bump. C IV 5808 
A is commonly seen in the red bump dKunth & Schildl 1 1986t 
IDinerstein & Shields|[T986b . 

The number of WR stars can be estimated from the lumi- 
nosity of the blue and red bumps divided by the luminosity of 
a sing le WR star in a certain subtype (see, e.g., IGuseva et al.l 
2000). In the following three subsections, we firstly attempt to 
identify the WR features in our GRB hosts, and then to estimate 



3.1. Identifying the WR features 

To be able to identify the WR features reliably, we select the 
high quality spectra of GRB hosts. The criterion for this selec- 
tion is assumed to be S /N > 10 in the continuum on both sides 
of the blue bump to ensure that the WR features detected are real 
and that a final selected sample is as large as possible. For other 
cases, there are very few chances to identify the WR features, 
even if they do exist in the spectra. This final selected sample 
contains 5 objects: GRB 980703, GRB 020405, GRB 020903, 
GRB 031203, and GRB 060218. 

We identify the WR features of these 5 GRB hosts based on 
the appearance of the blue bump and WR e mission lines by fit- 
ting th em. This fitting method is also used bv lBrinchmann et aU 
(2008). The fitting is performed using the SPECFIT routine in 
IRAF. We fit the whole blue bump with a single broad Gaussian 
profile. To fit the WR lines N V 4605, 4620 A, N III 4640 A, 
C III/C IV 4650, 4658 A, and He II 4686 A, when they exist, we 
adopt multiple Gaussian profiles to distinguish the broad com- 
ponent, which is mainly from WR stars, from the narrow com- 
ponent, which is mainly from nebulae. Nebular lines, blended 
within the bump, are also fitted. The FWHMs of narrow com- 
ponents of WR lines are set to be consistent with those of nebular 
lines. When carrying out the fit, the best-fit continuum is chosen 
very carefully, because the WR features depend sensitively on 
the continuum estimates. We limit the overall wavelength shift 
of the blue and red features to be \AA\ < 3A and the width of 
the broad component of WR lines to be ~ 3000 km s -1 FWHM , 
which are reasonable values found by (Brinchm ann et alj |2008). 
The fluxes of WR bumps and WR lines are measured from the 
fitted spectra. The fits are shown in Fig. [2] A significant He I 
5875 A emission line, evidence of young stars such as WR stars, 
is detected in most of objects. 



3.2. Identifying the subtypes of WR stars and estimating 
their numbers 

The WR features in galaxies originate mainly in two types 
of WR stars (WN and WC). In the WN types, the emis- 
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GRB980703 



Blue bumr 
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T 1 .0 

o 
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3 0.0 



4400 4600 4800 

Wavelength (A) 



GRB020405 
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Wavelength (A) 



, i , 1 1 i i 
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Fig. 2. The spectra of GRB hosts are plotted in the left panels. The continuums are marked using red lines, which are the results of 
the Starlight fit using BC03 templates. The blue boxes indicate the WR blue and red bumps. The WR emission lines are marked in 
red dot-dashed lines. In the right panels, the fits of WR bumps are plotted in red profiles. The fits of individual WR line and nebular 
line are plotted in green and blue profiles respectively (for more details, please see the online color version.) 



sion lines from helium and nitrogen ions are often seen in 
the spectra. The emission lines from helium, carbon, and 
oxygen can be considered as characteristic features of WC 
types. The relative strengths of these emission lines deter- 
mine the earl y (E) and late (L) subtypes of WN stars (WNE, 
WNL ~) (e.g.. Ivan der Hucht et alJ 11981b IConti. Leep & Perrvl 
1983HVacca & Contilll992h, and WC stars (WCE WCL) (e.g., 
Torres. Conti & Massevl ll986: Vacc a & Contill 19921) . In the blue 
bump, the broad WR emission lines, such as N V 4605, 4620 A, 



N III 4634, 4640 A, C III/C IV 4650, 4658 A, and He II 4686 A 
are mainly produced by WNL and WCE stars. In the red bump, 
C IV 5808 A is emitted by WCE. Normally, WNE stars can- 
not be distinguished from other WR stars, since they can emit 
all these lines. However, their contribution to the bump can be 
ignored because of their lower luminosity and shorter lifetime 
than WNL stars. WCL stars can be identified on the basis of 
their C III 4650, 5696 A emission lines. Another type of WR 
star, WO, can be identified using oxygen lines. Subtypes of WO 
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stars can be cl assified according to the relative strengths of thei r 
oxygen lines (iBarlow & Hummed fl98l IVacca & Contilll992l) . 
In our sample, no object exhibits the features of WO stars. 

The number of WR stars in a certain subtype can be esti- 
mated from the total luminosity of the blue and/or red bumps 
divide d by the luminosity of a single WR star dGuseva et al.l 
2000). To determine accurately the luminosities of the bumps, 
the contributions from nebular lines should be subtracted from 
the bumps. 

The number of WCE stars can be derived directly from 
the luminosity of the red bump, since that luminosity is only 
produced by WCE stars. WCE stars are usually represented 
by WC4 stars. The average C IV 5808 A line lumi nosity 
of LMC WC stars is 3.2 x 10 36 erg s" 1 (ISmithet al.l[l99oi: 
ICrowfher & Hadfieldll2006h . 

To determine the number of WNL stars, we consider the lu- 
minosity of the blue bump from which we subtract the contribu- 
tion of WCE stars. The contribution of WCE in the blue bump 
can be estimated from the ratio of the luminosity of WC4 stars 
in the red bump to that in the blue bump, 



^WC4 



(.14658) 



^WC4 



(i5808) 



(1) 



Several values of coefficient k were given by p revious works. We 
adopte d the value of k = 1.71 + 0.53 given bv lSchaerer & Vaccal 
(Il998h . WNL stars are usually represented by WN7 stars. The 
luminosity of a single WN7 in the blue bump, which depends on 
metallicity, is 2.0 x 10 36 ergs" 1 forZ < Z Q , or 2.6 x 10 36 erg s" 1 
for Z > Z (ISchaerer & Vaccal il998). The first is adopted in this 
work, since our samples have lower metallicity than solar (see 
Sect. 4). 

The number of WCL stars can be estimated from the lumi- 
nosity of C III 5696 A. For a WCL single star, the value of lumi- 
nosity of C III 5696 A adopted in this work i s 8.1x 10 36 erg s _1 
(ISchaerer & Vaccalll998l: iGuseva et al.ll2000h . However, the lu- 
minosity of this weak line in a single star is still not well known. 

The ratio of relative WR/O star number can be estimated 
by applyi ng two methods to optical spectra. The method devel- 
oped bv lArnault et al.l (|1989[) uses the flux in the entire WR blue 
bump, while another ( Vacca & Contil[T992T) uses only the flux of 
the broad He II 4686 A line. We use the first method in this work 
because it is more suitable for low resolution spectra, for which 
the broad He II 4686 A line cannot be separated clearly from the 
bump. The formula is 



log[ 



WR 



(WR + O) 



] = (-0.11+ 0.02) + (0.85 + 0.02) log ( Lbb " bump ). 

Lh/3 



(2) 

For each GRB host, both the fluxes and luminosities of the WR 
bumps and other weak WR features are given in Table[5] Fluxes 
given here are not corrected for interstellar extinction, while lu- 
minosities are absolute ones, which are transformed from these 
fluxes after being corrected for extinction. 

Several sources contribute to the uncertainties in estimating 
the numbers of WR stars. The luminosity of a single WR star 
adopted is based on the average of LMC and Milky Way WR 
stars, which could be higher than our host galaxy sample, be- 
cause of the difference in metallicity. The difference in metallic- 
ity could also affect the value of k. Besides these two major ef- 
fects, there are three other sources: the observational uncertainty, 
the errors in the multi-component profile fit, and the contamina- 
tion by nebular lines to the blue bump. Together, they produce 
the measured an uncertainty, assumed to be 30 - 60 %, in the 
number of WR stars and the WR/O star number ratio. 



3.3. The results 

We present the WR features and the number of WR stars in 
each of our GRB hosts in detail as follows, and as also given 
in Table 

GRB 980703 - The S/N in continuum on both sides of 
the blue bump (hereafter S IN*) is about 12. We do not detect 
the convincing WR bump in the spectrum. However, two WR 
lines, C IV 4658 A and He II 4686 A are identified, which may 
be emitted by WNL an d WCE stars (ISchaerer & Va cca 1998; 
iBrinchmann et alj|200 8). The red bump is not included within 
the spectral range of the spectrum. Therefore, we cannot verify 
the WR star subtypes. 

GRB 020405 - The S/N* is 10. A blue bump can be seen 
clearly in the spectrum. Some broad WR lines, N II 4640 A, 
C III 4650 A, C IV 4658 A, and He II 4686 A are identified in the 
bump, which implies the presence of WNL and WCE. The weak 
C III 4650 A indicates the presence of WCL stars. The WR sub- 
types and the number of WR stars cannot be verified since there 
is no red part to the spectrum. 

GRB 020903 - The S /N* is about 16. A strong blue bump is 
detected in the spectrum. Broad lines N V 4605, 4620 A, C IV 
4658 A, an d C IV 5808 A indicate the existence of WNL and 
WCE stars (IBrinchmann et al.ll2008h . The number of WCE stars 
is estimated by using the luminosity of the red bump, which can 
be interpreted as being produced by 1070 + 210 WCE stars. 
Adopting the value of k = 1.71 + 0.53, the contribution of these 
WCE stars in the blue bump is 5.49 + 1.81 x 10 39 erg s _1 . 
After subtracting the contribution of WCE stars, the luminos- 
ity of WNL stars in the blue bump is 4.23 + 1.46 x 10 39 erg s _1 , 
which corresponds to 21 17 + 721 WNL stars. 

GRB 031203 - The S/N* is about 20. The strong signature 
of the blue bump is shown in the spectrum. The broad WR lines 
C IV 4658 A, He II 4686 A are detected in the blue bump. In the 
red part, we do not detect the characteristic line of WCE, C IV 
5808, which implies that the number of WCE in this galaxy can 
be ignored. Besides that, the other two WR lines, N III 4905 
A, N II 5720-40 are also detected. All of these WR features are 
indicative of WNL stars dSmifh et al.lll996b IGuseva et al.l i2000: 
IBrinchmann et all 12008b . The number of WNL stars is 1585 + 
540 as estimated from the luminosity of the blue bump. 

GRB 060218 - The S/N* is higher than 40. Many WR 
lines are identified, including N III 4512 A, N V 4620 A, N II 
4640 A, He II 4686 A, He I/N II 5047 A, and N II 5720-40 
A. All these l i nes exhibit the character i stics of WNL stars 
(iMassevI l200H ICrowther & Smith! fl997t IGuseva et al.l l2000l: 
IBrinchmann etai1l2008l) . Two WR lines, C III 4650 A and C II 
5696 A, are also identified, which are characteristic of WCL 
stars (IBrinchmann et al.1l2008h . The presence of WCE stars can- 
not be confirmed, since no C IV 5808 line is detected. Using the 
luminosity of the blue bump, we estimate the number of WNL 
stars to be 324 + 120. The luminosity of C III 5696 A corre- 
sponds to 31 + 10 WCL stars. 

The relative WR/O star number ratios are estimated using 
the luminosities of blue bumps and H/3. For the hosts of GRB 
020405, 020903, 031203, and 060218, the ratios are 0.04+0.03, 
0.21 +0.06, 0.01 +0.01 , and 0. 1 1 +0.04, respectively. 

In summary, obvious WR features are detected in 4 out of 5 
spectra of long-duration GRB hosts, except for GRB 980703, 
which shows no convincing WR bump but two WR lines. 
Moreover, the numbers of WR stars for 3 GRB hosts are in 
the range from 300 to 6000, which is co nsistent with the typ- 
ical number of WR stars (~ 100 - 10 s . IVacca & Contil ll992: 
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Ion 


980703 


020405 


020903 


031203 


060218 




F L 


F L 


F L 


F L 


F L 


blue bump (A4650)" 
red bump (,15008) 
CIIL15696 




0.18±0.09 654.39±326.21 


1.52+0.31 97.22+42.90 
0.50±0.10 32.10±7.30 


10.65±6.06 114.87±65.37 


8.08+3.25 6.49+2.41 
3.76+1.10 2.48+0.76 


WNL 
WCE 
WCL 
WR/O 




0.04+0.03 


2117+721 
1070±210 

0.21±0.06 


5743±3268 
0.01+0.01 


324±120 

31±10 
0.11+0.04 



Notes: WR emission line fluxes (F, corrected for Galactic extinction) are in units of 10 17 erg s ' cm 2 , Luminosities (L) are in units of 
10 38 erg s _1 . (a) Fluxes and luminosities of nebular lines were subtracted. 



iGuseva et alJ l2QO0h in WR galaxies. High WR/O star number 
ratios are also estimated in these 4 GRB hosts Therefore, the 
presence of WR stars in these GRB hosts is verified. However, 
because of the small sample studied here, the conclusion that 
WR stars exist in all the long-duration GRB hosts still cannot be 
confidently made. The results of this work in all cases support 
the scenario that WR stars are the progenitors of long GRBs, 
and are consistent with th e previous studies described below. 

lHammer et al.1 (120061) first discovered WR stars in three 
nearby long-duration GRB host galaxies (GRB 980425, GRB 
020903, GRB 031203) using deep spectroscopic observations. 
Since then, little ha s been done to detect WR features in long- 
duration GRB hosts. Margutti etalJ(l2007l) claimed that the iden- 
tification of WR emission lines in GRB 031203 is uncertain. 
However, their Fig. 3 does show the strong WR blue bump. 
The N III 4640 A and C IV 4658 A blended broad lines can 
also b e identified in the bum p, which is not attempted in their 
work. IWiersema et alJ (|2007) did not detect the obvious He II 
4686 A and WR bump in GRB 060218. However, we note that 
the UVES spectra that they used were taken at a time close to 
the burst, which means that the significant contamination by the 
associated SN 2006aj cannot be ignored in the WR feature iden- 
tification. In their work, an upper limit to the ratio WR/(WR + 
O) for the galaxy is estimated to be < 0.4, which is consistent 
with the result (WR/O = 0.1 1+0.04) estimated in this work. 
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Fig. 3. The diagnostic diagrams of GRB host galaxies. In the left 
panel , [O III] fH0 vs. [N II]/Ha relations is shown. The relation 
from lKewlev et al.l(l2001l) is drawn in the d ot-dashed curve. The 
solid curve is fro mlKauffmann et"aD(l2003h . [O III] /H/J vs. [S II] 
/Ha relation from Kauff mann et al.l d2003l) . is shown in the right 
panel. 

that the GRB host galaxies are star-forming galaxies and AGN 
contamination in GRB host galaxies can be ignored. 



4. The physical properties of GRB host galaxies 

Metallicity is a fundamental parameter for probing the proper- 
ties of GRB progenitors and the environment of GRB regions. 
Strong emission lines in spectra of GRB host galaxies allow us to 
estimate metallicities. The luminosity (stellar mass)-metallicity 
(L - Z, M, - Z) relation of galaxies is a fundamental relation for 
indicating their evolutionary status and star-formation histories. 
We estimate the metallicities of these GRB hosts, and then study 
them in terms of the L - Z and M» - Z relations. 

4.1. AGN contamination 

The traditiona l Baldwin-Phillips-Terlevich (BPT) diagram 

dBaldwin et al.l Il98ll) can diagnose the origin of the narrow 
emission lines of emission-line galaxies. Fig. [3] shows the lo- 
cations of GRB host galaxies in the BPT diagra m. The [O III] 
/H/j vs . [N III/Hor diagnostic lines a re taken from lKewlev et al] 
d200l and iKauffmann etafl J20031 T he [O III] /H/J vs. [S II] 
/Ha relation from lKewlev et al.l (12001 1) is shown. This illustrates 



4.2. Metallicity 

In terms of oxygen abundance, the metallicities of GRB hosts 
can be estimated using the T e method and R23 methods. The 
T e method based on electron temperature is a "direct" way. The 
"indirect" method is the so-called "strong-line" method, such as 
the /?23 method based on the empirical relationship between O/H 
and 7?23 ■ 

It is well known that the T e method is effective for metal- 
poor galaxies because the characteristic line [O III]4363 is only 
observable in this case. Out of 8 galaxies, 5 GRB hosts (GRB 
990712, 020903, 030329, 031203, and 060218) clearly exhibit 
an [O III]4363 line, which implies that they have low metal- 
licities and lie on the lower-metallicity branch of the R23 di- 
agnostic diagram. [N II]/[0 II] also provides a r eliable means 
of loca ting a galaxy on the metallicity branch. As iKewlev et al.l 
(2008) suggested, the galaxies with log([N II]6583/[0 IIJ3727) 
<-1.2 lie on the lower metallicity branch. All of these 5 
objects and the GRB region of GRB 060505 have values 
of log([NII]6583/[OII]3727) <-1.2, while the entire GRB 
060505 host galaxy lies on the upper branch, because of the 
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value of its log([N II]6583/[0 II]3727) >-l.2. The values of 
log([NII]6583/[OII]3727) are given in Table For another 2 
GRB hosts, 980703, 020405, we cannot verify the metallicity 
branch which on they lie because their [NII]6583 lines are not 
detected. Moreover, we use log([N II]6583/Ha) to verify the 
metallicit y branches (<- 1.3 for lower branch; for upper 

branch. iKewlev et"afll2008l) for the sample. These three indica- 
tors provide the consistent results for these galaxies. 

For the 5 host galaxies that exhibit the [O III]4363 emission 
line, we use the T e method to calculate their oxygen abundances. 
Electron temperature and electron density are estimated using 
the TEMDEN task of IRAF, which is based on a 5-l evel atom 
program and was described by IShaw & Dufo ur ( 1994). The flux 
ratios of [O III] 7(4959 + 5007)//(4363) are adopted to calcu- 
late electron temperature when available. Using the flux ratio 
of [S II] 7(6716)//(6731), the electron density can be derived, 
which is needed in the electron temperature and oxygen abun- 
dance calculation. 

When carrying out the electron density calculation, the [S II] 
7(6716)//(6731) doublet is not available or too noisy for some 
objects in our sample. Therefore, we tested the sensitivities of 
both electron temperature and oxygen abundance to electron 
density, finding that neither is very sensitive to electron den- 
sity. This allows us to estimate oxygen abundance for a wide 
range of electron density (20 - 200 cm 4 ) without significant bias 
(the discrepancy is less than 5 %). To calculat e TVbased oxyge n 
abundances, we us e the methods des c ribed by Stasiriska (2005), 
llzotov et al.1 d2006l) . lYin et alj (120071) . lEiang et all (12006I) . 

The 7?23 method has b e en extensively di s cussed in the lit- 
erature dPagel er all 1 1 979b iMcGaughj 1 199 li iKobulnickv et all 
11999b iTremonti et al.1 120041: fan et al.N2007l: Liang et al.ll2007l) . 
To estima te their i?23-based metallicities, we adopt the for- 
mulae of Kobulni ckv et al.1 (1 19991) for both the upper and 
lower branches of /?23-(0/H) solutions. We use indicators of 
[N II]6583/[OII]3727 and [N II]6583/Ha to break the R 2i de- 
generacy. The estimated oxygen abundances and errors are listed 
in Table [6] The r c -based and R23 -based metallicity estimates 
have an acceptable discrepancy of 0.03-0.3dex. 

We now compare our metallicity estimates with those in the 
literature. There are four objects for which 7Vbased oxygen 
abundances ha ve been estima t ed in the literature. For the host of 
GRB 020903. iHammer et all (120061) measured a 12+log(0/H) rf 
of 7.97, which is very consistent with this work (8.01) within 
the errorbars; Savaglio09 gave 8.22, which is 0.21 d ex higher. 
For the host of GRB 030329, iLevesque etail (2009) obtained 
7.72 from Keck spectra, while we obtain 7.65. These two es- 
timates are extremely consistent within the er rorbars. For the 
host of GRB 03120 3, iProchaska et ail d2004 obt a ined 8 .10, 
iMargutti et ail d2007l) inferred 8.12, ILevesque et all d2009l) ob- 
tained 7.96, and Savaglio09 measured 8.02, while we obtain 
8.14. All of these estimates are con sistent within the error bars. 

For the host of GRB 060218. Iwiersema et all d2007l) m ea- 
sured a 12+log(0/H) value of 7.54. ILevesque et al.1 d2009l) ob- 
tained 7.62, Savaglio09 measured 7.29, and we inferred 7.88, 
which is closest to the second value. The discrepancy between 
these works could be caused by contamination by the associ- 
ated supernova. Savaglio09 gathe red t he very early-time spec - 
tra observed by [Pian et al.l d2006l) and ISollerman etail d2006l) . 
IWiersema et al.l d2007l) acquired their spectrum within 1 month 
of the burst, while our spectrum of the host was taken in Dec. 
2006, 10 months after the burst, and lLevesque et all (12009) took 
the spectrum even later (19 months after the burst). In this case, 
the contamination of the GRB 060218 host galaxy by the light 
of SN2006aj was significant and lasted for several months. This 



Table 8. Magnitudes and stellar masses of GRB host galaxies 



GRB 



M K 



logM,[M Q ] 



980703 
990712 
020405 
020903 
030329 
031203 
060218 
060505" 



-20.85 ±0.12 
-18.67 ±0.09 



-18.96: 
-16.22: 

-20.58 : 

-16.02: 
-19.24: 



0.10 
0.07 
0.05 
0.07 
0.02 



-24.41 : 
-21.85 : 



-22.27 : 
-17.23 : 
-22.20 : 



0.12 
0.10 



0.02 
0.34 
0.04 



11.14: 
10.14: 



0.39 
0.34 



10.13 ±0.26 
8.28 ±0.55 
10.18 ±0.24 



Notes: Absolute magnitud es M B and M K are corrected for galactic 
foreground extinction ( Schlegel et al. 1998). (a) The magnitudes 
and stellar mass of the entire galaxy of GRB 060505 host 



contribution strongly affects measurements of emission lines; 
hence, the errors in the metallicity estimates for the early spec- 
tra mi ght be large. Our result is close to that of ILevesque et al.l 
d2009l) . which is also derived from the late-time observation. 
Anyway, all of them show this host has low metallicity. 

We do not compare our R^i -based estimates with any other 
work, because different calibrations provide very different val- 
ues, and we discussed the discrepancy between our own T e and 
7?23-based abundances above. 



4.3. Luminosity vs. metallicity 

We obtain multi-band photometry from literature, which is 
shown in Table Q All apparen t magnitudes are co rrected for 
Galactic foreground extinction (Schle gei et al] [i~998) and inter- 
nal extinction. The k-correction is performed fo r the magnitudes 
toz = by using the kcorrect v4- 1 -4 programp| (jOke & Sandage 
119681; iHogg et alJl2002t iBlanto & Roweisll2007l) . The absolute 
B- and K-bwd magnitudes corrected for Galactic extinction are 
given in Table [8] We plot our long-duration GRB hosts on the 
L - Z diagram (Fig. 0), which is presented in terms of ab- 
solute B magnitude and the /?23-based metallicities. The mag- 
nitude and metallicity of the entire GRB 060505 host galaxy 
are presented in the L - Z diagram. For GRB 980703, the 
mean value between the lower and upper branch is adopted, 
when the 7^23 degeneracy cannot be broken. The L - Z rela- 
tions of various low and high redshi ft samples from t he lite rature 
(SDSS star-forming gal axies from Tremonti et ail d2004l) . UV- 
selected galaxies fro m IContini et al. d2002l). l arge magnitude- 
limited sample from Lamareille et al.l (| 2 04). emission-line- 
selected galaxies from Melbourne & Salzer (2002), irregular and 
spiral galaxies from Kobulnickv & Zaritsky] (1 19991). and irreg- 
ular g alaxies from Skill man et al.l d 19891) and iRicher & McCalll 
(fl995h are also shown in Fig. |4] In comparison with the L - Z 
relation, the GRB hosts show an obvious discrepancy from 
other samples. GRB hosts have lower metallicity values than 
other galaxies at given luminosities. For luminous GRB hosts, 
this discrepancy is e ven larger. This trend is also shown in 
ILevesque etail d2009l) . 



4.4. Stellar mass vs. gas phase metallicity 

The stellar mass of galaxy can be e stimated f r om th e stellar 
mass-to-light ratio and color following lBell et akl d2003l) . We use 
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Table 6. Metallicities of GRB host galaxies 
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980703 


990712 


020405 


020903 


030329 


031203 


060218 


060505 s 


060505'' 


N e [SII] 11 




20 ~ 200 




20 ~ 200 


20 ~ 200 


52.63 ± 43.20 


120.68 ± 32.22 






T c [OIII] b 




16434.9 ± 425.8 




13517.1 ± 335.2 


17911.7 ±325.0 


12902.2 ± 305.2 


15113.3 ± 335.5 






12 + log (O/H) 




















(Tef 




7.812 ±0.15 




8.014 ±0.15 


7.653 ±0.14 


8.142 ±0.07 


7.883 ± 0.09 






12 + log (O/H) 




















(K99J) d 


8.28 ±0.14 


8.18 ±0.12/ 


8.10 ±0.12 


8.11 ±0.12' 


8.06 ±0.10'' 


8.11 ±0.11^ 


8.16 ±0.14' 


8.06 ±0.14' 


7.99 ± 0.09 


12 + log (O/H) 




















(K99.u) d 


8.44 ±0.15 


8.47 ±0.15 


8.53 ±0.15 


8.50 ±0.15 


8.56 ±0.16 


8.51 ±0.15 


8.50 ±0.16 


8.57 ±0.15 


8.66 ±0.16^ 


log([NII]6583 




















/[OII]3727) 




-2.36 




-1.47 


-2.12 


-1.23 


-1.45 


-1.36 


-0.52 


A V 


0.144 ±0.532 


0.551 ±0.025 


1.927 ±0.359 


0.146 ±0.086 


0.244 ± 0.046 


0.280 ± 0.040 


0.000 ±0.018 


0.526 ±0.081 


1.601 ±0.310 



Notes: (a) Electronic density estimated from the [S II] flux ratio: 7(6716)//(6731). (b) Temperature in K estimated from the following flux ratios: [O III] 
7(4959 + 5007)/7(4363 ) (c) Metallicity, 12+log(0/H), estimated using the effective temperature method, (d) Metallicity, 12+log(0/H), estimated following 
Kobulnicky et al. 1 1999). K99J (their Equation 8) is valid for the lower branch; K99_h (their Equation 9) is valid for the upper branch, (e) Extinction coefficient, 
Ay (in magnitude) is derived using the standard Balmer ratio of Ha and H/3. For GRB 980703 and 020405, Ay from the ratio of Hy and H/3 is adopted, (f) 
Adopted metallicity in this work. The R23 degeneracy are broken using the [N II]/Hq' and [N II]/[0 II] ratios, (g) The GRB site (h) The entire host galaxy 



Table 7. Photometry of GRB host galaxies 



GRB 


U 


B 


V 


R 


1 


J 


H 


K 


Refs 



980703 - 23.40±0.12 23.04±0.08 22.58 ± 0.06 21.95 ±0.25 20.87 ±0.11 20.27 ±0.19 19.62±0.12 1,2,3 

990712 23.12 + 0.05 23.36 ± 0.09 22.39 ±0.03 21.84 ±0.02 21.41 ±0.03 20.81 ±0.17 20.25 ±0.19 20.05 ±0.10 3 

020903 - 21.70 ±0.10 20.80 ±0.10 20.80 ±0.10 20.50 ±0.10 - - - 4,5 

030329 22.68 ±0.10 23.42 ±0.07 22.88 ±0.05 22.80 ±0.04 - 21.52 ±0.04 21.17 ±0.24 6 

031203 22.41 ±0.18 22.32 ± 0.05 20.53 ± 0.05 20.44 ± 0.02 19.40 ±0.04 18.28 ±0.02 17.78 ±0.02 16.54 ± 0.02 7 

060218 20.45 ±0.15 20.46 ± 0.07 20.19 ±0.04 19.86 ± 0.03 19.47 ±0.06 18.99 ±0.16 18.52 ±0.22 18.69 ±0.34 8,9,10 

060505 18.43 ±0.05 18.89 ± 0.02 18.27 ±0.02 17.90 ± 0.02 17.51 ±0.02 17.29 ± 0.08" - 15.85 ± 0.04 11 

F555W a F102W - 

020405 - - 22.63 ±0.05 21.84 ±0.05 21.29 ±0.05 - - - 12 



References: (1 ) IVreeswiik et al.) dl999h : ( 2) ISokolov et all d200lh : (3) IChristensen et alJ d2004h: (4) ISoderberg et al] (120041) : ( 5) iBersier et al. 
( 20061) : (6) lGorosabel et al.1 d2005h ; ( 7)lMargutti et alj d2007hr (8) ISollerman et all d2006l) ~(9) ICastro Ceron et alj d2008h : (lO) lKocevski et al. 



(2007D; (ll) lThQne et alJ d2008): (12) Wai nwright et all d2007 
Notes: The magnitudes are not corrected for Galactic extinction, (a) F555W (Central wavelength: 5407 A; Band width: 1236 A); F702W (Central 
wavelength: 6895 A; Band width: 1389 A); F814W (Central wavelength: 7940 A; Band width: 1531 A) 



the (B - V) color and Mk to estimate the stellar masses of GRB 
host galaxies. The formula is 

log A = -0A(M K - 3.28) + [a K + b K (B -V) + 0.15], (3) 

where Mk is the /T-band absolute magnitude and (B - V) is the 
re st-frame color. Th e coefficients ok and come from Table 7 
of lBell et al] (120031) . The stellar masses are listed in the last col- 
umn of Table [8] We compare our results with previous work on 
GRB hosts dCharv et al.l 120021: ICastro Ceron et aLll2006l 120081 
ISavaglioet al.ll2009h . Our results are more consistent with the 
stellar masses derived by ICastro Ceron et al] ((2008) within the 
errorbars. However, our results are systematically higher than 
theirs by ~ 0.3 dex. This discrepancy between the two datasets 
could be caused by the use of different M*/Lk ratios. The stel- 
lar m asses given by Savaglio09 are the lowest of all previous 
work. lCastro Ceron et al.l d2008l) explain that the reasons for this 
discrepancy could be the different M»/L# ratios applied and the 
underestimated dust extinction in Savaglio09, which is also con- 
firmed by our dust extinction (see TableHJl. 

We plot our long-duration GRB hosts on the stellar mass- 
metallicity (M* - Z) diagram (Fig. |5). The T e - and /?23-based 



metallicities are adopted and plotted in the left and right pan- 
els, respectively. In the right panel, the open circle and black 
filled star represent the hosts of GRB 060505 (not the GRB re- 
gion) and GRB 980703, respectively. The M, - Z relation of 
lo cal normal star-fo rming galaxies selected from SDSS derived 
bv lLiang et al] d2007l) is given in Fig.[5](the red dashed line de- 
rived from the T e method in the left panel; the red dotted line 
derived from the R23 method in the right panel). The T e - and 
/?23-based metallicities of long-duration GRB host galaxies are 
all obviously lower than those of local star-forming galaxies with 
comparable stellar masses. 

5. Discussion and conclusions 

We have attempted to identify the WR stars in GRB host galax- 
ies. According to the collapsar model, WR stars are consid- 
ered as the most favored candidates of the progenitor of long- 
duration GRBs, therefore the presence of WR stars in long- 
duration GRB hosts can directly validate this model. We have 
investigated a sample of 8 GRB host galaxies from archival data 
from VLT/FORS2. Out of those 8 GRB hosts, 5 galaxies hav- 
ing spectra with S /N > 10 were for our sample to detect the 
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Fig. 5. In the left panel, the r e -based (the filled circles) oxygen abundances of long-duration GRB h osts are given. Th e red dashed 
line refers to the M* - Z relation of local star-forming galaxies from SDSS derived from T e method dLiang et al.ll2007l). In the righ t 
panel, we plot the 7?23-based oxygen abundances of long-duration GRB hosts derived using formulae from Kobulnickv et al. ( 1999) 
(the blue filled circles - oxygen abundances that break the degeneracy; the black star - oxygen abundances from mean value between 
the lower and upper branch). The open circle represents the entire host of GRB 060505. The re d dotted line refers to the M, - Z 
relation of local star-forming galaxies from SDSS derived from R23 method dLiang et alj |2007). See the online color version for 
more details. 



WR features. The presence of WR features was verified in this 
entire sample. Out of those 5 GRB hosts, 4 certainly have WR 
features detected, which are GRB 020903, GRB 031203, GRB 
020405, and GRB 060218. For the remaining one, GRB 980703, 
the detection is marginal. The subtypes and numbers of WR stars 
in those 5 GRB hosts were also derived. Those results strongly 
support the collapsar model by illustrating the link between WR 
stars and GRBs. 

Additional aim was to investigate the physical properties of 
long-duration GRB host galaxies, such as metallicity, luminos- 
ity, and stellar mass. Comparing with literature, we used the 
most consistent measurements and calculation methods for all 
8 objects to ensure systematic coherence and comparability. We 
found that the long-duration GRB hosts show obvious disagree- 
ment with the L-Z and M» - Z relations derived for low redshift 
galaxies in several samples. The luminous and massive long- 
duration GRB hosts have lower metallicities than other galaxies 
in these samples. 

The WR/O star ratio in WR galaxies is found to decrease 
as a fu nction of decreasing metallicity bv ICrowther & Had field 
(2006). The WR/O star ratio and metallicity in GRB host galax- 
ies do not show a clear trend. However, compared to WR galax- 
ies dGusev a et al. 2000), the observed WR/O ratios are higher in 
the GRB hosts studied here, while metallicities in GRB hosts are 
obviously lower than in the numerous population of low redshift 
galaxies. This suggests that the hosted regions of GRBs are con- 
sistent with the first stage of star formation in a relatively pristine 
medium. Furthermore, the ratio of WC/WN stars is also found 
to decrease with metallicity. The relation betw een the ratio o f 
WC/WN stars and metallicity was derived by Massey (2003). 
The hosts of GRB 020903, 060218 deviate from the relation, in 
having higher WC/WN ratios of ~ 0.5 and 0.1 than Local Group 
galaxies (see their Fig. 1 1 , the upper panel). The collapsar model 
suggests that the WC stars are more likely to be the progenitors 



of long-duration GRBs than the WN stars. Therefore, the ob- 
served high WC/WN ratio could also be evidence to support this 
model, although the sample is small. 

One of our galaxies, the host of GRB 060505, is well re- 
solved in the ground-based observation. However, the type of 
this burst remains a topic of debate. GRB 060505 has a burst 
duration of ~ 4 s, but lacks evidence of an accompanying super- 
nova. It is classified as a long-duration G RBs that have no as - 
sociated s upernovae (e.g., GRB 0606 14) bv lThone et al] d2008l) , 
However, iLevesque & Kewlevl d2007l) suggest that the environ- 
ment of GRB 060505 is more consistent with the host environ- 
ments of short-duration GRBs. We have investigated the metal- 
licity of the GRB site and the entire host galaxy separately. 
We have found a relatively low metallicity in the GRB region 
and a highe r one i n the entire host, which are consistent with 
iThone etafl d2008l) . Unfortunately, the S /N of its spectrum is 
not high enough to detect WR star features. However, a very 
young (~ 6 Myr) s tellar population in the GRB site is found by 
IThone et alJ (|2008) using stellar population modeling. This low 
age corresponds to the lifetime of a 32 M G star. Evidence from 
the properties of the GRB region suggests that the GRB 060505 
originated in a long-duration core-collapse progenitor. 

Most of the GRB host galaxies in our sample are too dis- 
tant to allow spatially resolved analysis with ground-based spec- 
troscopic observations. Therefore, the only information that we 
can derive concerns about the global properties of galaxies. The 
properties of the actual explosion sites are still unclear. For 
further study of GRB progenitors and properties of GRB host 
galaxies, more high-resolution images and deep spectroscopy 
are needed. X-shooter at the VLT will be the ideal instrument 
to investigate in greater detail the chemical and stellar popula- 
tion properties of GRBs. 
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Fig. 4. The L - Z relation of long-duration GRB hosts. The 
/?23 -based metallicities are adopted. The blue filled circles repre- 
sent oxygen abundances which break the degeneracy of R23; The 
black star represents the mean value between the lower and up- 
per branch of GRB 980703. The open circle represents the entire 
host of GRB 060505. Luminosity-metallicity relation for SDSS 
galaxies and various galaxy samples are drawn from the litera- 
ture (see legend). See the online color version for more details. 
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